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ABSTRACT 

The characteristic size of early-type galaxies (ETGs) of given stellar mass is observed to increase 
significantly with cosmic time, from redshift z>2 to the present. A popular explanation for this size 
evolution is that ETGs grow through dissipationless ("dry") mergers, thus becoming less compact. 
Combining N-body simulations with up-to-date scaling relations of local ETGs, we show that such 
an explanation is problematic, because dry mergers do not decrease the galaxy stellar-mass surface- 
density enough to explain the observed size evolution, and also introduce substantial scatter in the 
scaling relations. Based on our set of simulations, we estimate that major and minor dry mergers 
increase half-light radius and projected velocity dispersion with stellar mass as R c oc M^ 09±0 - 29 and 
er C 2 oc JVfO° 7 ± 011 j respectively. This implies that: 1) if the high-z ETGs are indeed as dense as 
estimated, they cannot evolve into present-day ETGs via dry mergers; 2) present-day ETGs cannot 
have assembled more than ~ 45% of their stellar mass via dry mergers. Alternatively, dry mergers 
could be reconciled with the observations if there was extreme fine tuning between merger history and 
galaxy properties, at variance with our assumptions. Full cosmological simulations will be needed to 
evaluate whether this fine-tuned solution is acceptable. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: kinematics 
and dynamics — galaxies: structure — galaxies: evolution 
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1. INTRODUCTION 

In the hierarchical model of structure formation, 
galaxy mergers are expected to play a central role in 
the assembly and evolution of galaxies. If present-day 
early-type galaxies (ETGs) have experienced mergers in 
relatively recent times (z < 1 — 2), most of these mergers 
must have been dissipationless o r "dry" , because the old 
stellar populations of ETGs (e.g. lThomas "et~aT][2005l) arc 
inconsistent with substantial recent star formation. 

Present-day ETGs satisfy several empirical correla- 
tions between half-light radius R c , central stellar ve- 
locity dispersion cto, bulge luminosity L, stellar mass 
M* , black- hole mass A/bh a nd gravitational lensing mass 
Mi e ns, such as the L-oq dFaber fc Jackson! I1976D, L - 
R e (iKormendvl 119771). Mr H -£ (iMagorrian et al.l I1 998L 
Mrh ^ dFerrarese fc Merrittl 120001 : iGebhardt et all 
l2000r i. Mrh-M, (iMarconi fc HTrntT 120031). M \ e ^-R 
and Mi ens -(To relations (jNipoti. Treu. fc Bolton! 2009, 
hereafter NTB09), the Fundame ntal Plane (he r eafter 
FP, relating L, R e an d a ; iDressler et al.l Il987t 
iDiorgovskv fc Davis! Il987h and the lensing Mass Plan e 
(relating Mi ens , <J and R e ; iBolton et all I2007L [2008b] l. 
The existence of these scaling laws has been used to 
establish that present-day ETGs did not form by sev- 
eral dry mergers of lower mass progenitors obeying the 
same scaling laws, because as an effect of dry mergers the 
galaxy size increases rapidly with mass, while the veloc- 
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ity dispersion remains almos t constant, in contrast with 
the observed correlations dCiotti fc van Albadal 120011: 
Nipoti, Londril lo. fc Ciottil l2003t iBovlan-Kolchin et al.l 
2006; Ciotti. Lanzoni fc Volonterill2007l NTB09). 

Photometric and spectroscopic observations of high- 
rcdshift (z > 1 — 2) ETGs suggest that these objects 
may be remark ably more compact t h an their local coun 
terparts (e.g. iStiayelli et al.l {l99l iDaddi et alJ ] 2Qi 
Truiillo et al.ll2006HZirm et al.ll2007HCimatti et al.l l2& 
van der Wei et all l2008t Ivan Dokkum et all 120081 here- 



after vD08; Saracco, Longhetti, & Andreon 2009). These 
findings contributed to renew interest in the proposal 
that dry merging might be a key process in the growth 
of ETGs, because it is one of the few kn own mechanisms 
able to make galaxies less compact fe.g.lKhochfar fc Silk 



2006 ; lHopkins et"aL 2009a: Naab . Johansson, fc Ostriker 
2009t Ivan der Wei et al.ll2009D . However, although dry 
merging works qualitatively in the right direction, it is 
not clear whether it works quantitatively. 

In this Letter we address the question of the viability 
of dry merging as a process driving the structural evo- 
lution of ETGs, by testing whether it is able to produce 
the local scaling laws with their relatively small intrinsic 
scatters. For this purpose, we use the set of numerical 
simulations of dissipationless minor and major mergers 
presented in NTB09 and test them against the stellar- 
mass scaling relations of the S loan Lenses ACS Survey 
(SLACS) sample of lens ETGs (lAuger et al.H2009D . This 
strategy allows us to construct a direct link between the 
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observed properties of high-redshift {z > 1— 2) galaxies — 
for which estimates of the stellar masses are available — 
and realistic simulated galaxies, in which the relative 
distribution of luminous and dark matter has been cali- 
brated on the lensing-mass scalin g relation s (NTB09). 

Throughout the paper we adopt lSalpeterl (| 19551 ) Initial 
Mass Function (IMF), which is consistent wit h the gravi- 
tation al lensing constraints for SLACS ETGs (jTreu et al.l 
2009b). Observational data are converted to this calibra- 
tion when necessary. The main results of this paper do 
not depend on this choice, since they involve stellar-mass 
ratios and the IMF normalization effectively factors out. 

2. OBSERVED STELLAR-MASS SCALING RELATIONS 

In the present work we consider the stellar-mass-size, 
stellar- mass-velocit y-dispersion a nd st ellar-mass FP re- 
lations derived by [Auger et all (|2009f ) for the SLACS 
sample of lens ETGs, assuming Salpeter IMF stellar 
masses and using the velocity dispersions and rest-frame 
V-band half-light radii derived from Sloan Digital Sky 
Survey (SDSS) spectroscopy and Hubble Space Tele- 
scope photometry. The correlations are obtained us- 
ing galaxies with stellar masses generally in the range 
11 < logM»/M <12. 

The stellar-mass-size relationship is given by 



l0g (lpl:) = (°- 730±0 - 047 ) 



x log 



10 n M Q 



0.391 ±0.029, 



(1) 



with intrinsic scatter in log R e at fixed mass 0.071. The 
stellar-mass-velocity-dispersion relation is 



logf , Ue2 A = (0.196 ±0.033) 
V kms _i / 



x lop 



2.281 ±0.021, 



(2) 



where cr e 2 is the projected velocity dispersion within an 
aperture radius R c /2. This has an intrinsic scatter in 
log<7 e 2 at fixed mass of 0.055. Finally we consider the 
FP-like correlation 



logc* = (0.099 ±0.067) 



x log 



10 n M Q 



0.959 ±0.042, 



(3) 



where c» = a^ 2 R c /^GM^ is a dimensionless structure 
parameter. The intrinsic scatter in logc* at fixed M* is 
0.077. 

We recall that SLACS lenses are indistinguishable from 
control samples of SDSS galaxies with the same stel- 
lar velocity dispersion and size, in terms of luminos- 
ity and distribution of environments ( Bolton et al.ll2006l : 
iTreu et al.|[2006l: iBolton et al.ll2008at iTreu et al.l l2009a) . 
In addition, as discussed in lAuger et al.l ( 20091 ). the above 
scaling relations are consistent with other scaling rela- 
tions derived for larger sam ples of SDSS galaxies wi th a 
similar range of masses (e.g. iHvde fc Bern ardi 200jl): in 
particular the intrinsic scatter of the M*-i? relation is 
consistent with the observed sc a tter o f the widely used 
correlation found by iShen et all (|2003f ) . 
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Fig. 1. — Structural parameter c* = ir^ 2 R a /2GM f ) (top panel), 
effective radius (central panel) and projected velocity dispersion 
(bottom panel) as functions of the stellar mass. The solid and 
dashed lines represent the local observed correlations (equations [1} 
with the associated intrinsic scatters. The big hexagon repre- 
sents the median M» and R a for a sample of nine z ~ 2.3 ETGs 
(vD08). The small hexagons with error bars represent the ETG 
1255-0 (z = 2.186; vDKF09). The other points and error bars 
represent the average values and the l-cr scatters (due to projec- 
tion effects) for the simulated galaxies of the merger hierarchies 
indicated in the top panel (/* is the stellar mass fraction, L is 
the orbital angular momentum modulus and N B is the number of 
satellites per step). 

3. EFFECT OF DRY MERGERS ON COMPACT HIGH-Z 
EARLY-TYPE GALAXIES 

Here we combine the scaling laws described above, ob- 
servations of high-z ETGs and the results of the N-body 
simulations presented in NTB09 to address the ques- 
tion of whether dry mergers can be responsible for the 
observed size evolution of ETGs. The set of simula- 
tions (see NTB09 for details) consists of both major and 
minor dry-merger hierarchies, i.e. sequences of two or 
three steps in each of which the galaxy mass increases 
by a factor of ~ 2 as a consequence of either a binary 
equal-mass merger or the accretion of 7V S smaller satel- 
lites. ^The_Sj3ed_gjaL^ are spherical, 
with lNavarro. Frenk &; White! ((1996) dark-matter distri- 
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Fig. 2. — The shaded areas represent the distributions in the M*- 
c*, M t -R e and M t -(J e 2 planes of allowed dry-merging progenitors 
of local ETGs with stellar mass Ml (vertical bars). The upper 
and lower boundaries of the shaded areas correspond to 1-tr of the 
distributions, while the solid lines within the shaded areas show 



their average values. 



R' e and a'„ are the average values of c», 



Re, and a e 2 for local ETGs with mass M' t 
correlations are plotted as in Fig. [T] 



The local observed 



bution, IDehnenl (|1993f ) stellar distribution (with central 
logarithmic slope 7 = 1.5), and no gas. Altogether we 
consider eight major-merger hierarchies (differing in the 
orbital energy and angular momentum of the galaxy en- 
counters, and in the stellar mass fraction /» of the seed 
galaxies) and four minor-merger hierarchies (differing in 
/* , N s and in the phase-space distribution of the satel- 
lites' centers of mass). Each galaxy collision is followed 
up to the virialization of the merger remnant. In the sec- 
ond and third steps of the hierarchies the remnants of the 
previous steps are used to realize the initial conditions. 

The seed galaxies (and consequently all the end- 
products of each merging hierarchy) can be rescaled by 
fixing arbitrarily the simulation mass and length scales 
M* and r* (see NTB09). Though the relative effect of 
dry mergers on the galaxy structure and kinematics is 
independent of M* and f*, it is useful to fix these quan- 
tities, without loss of generality, to make the comparison 



with observations straightforward. For our purpose, it is 
natural to consider seed galaxies with stellar masses and 
sizes similar to those of high-z ETGs. Here we take as 
reference the sample of nine z ~ 2.3 ETGs studied by 
vD08, having median stel lar mass (M f ) ~ 2 .7 x 10 n M Q 
(assuming Salpeter IMF; iKriek et al.ll200ct l and median 
effective radius (R c ) = 0.9 kpc. These galaxies have a 
factor of ~ 5 smaller size than local ETGs of similar 
mass, as apparent from the central panel of Fig.[TJ where 
their median location in the M*-R c plane can be com- 
pared with the stellar-mass-size correlation (fTJ observed 
locally. 

We then fix M* and f* such that all our seed galaxy 
models have total stellar mass M* = 2.7 x 10 11 M© 
and R c = 0.9 kpc (symbols inside the big hexagon 
in Fig. [TJ. As a consequence, the velocity unit is 
given, so our seed galaxy models have projected ve- 
locity dispersions a c i — 532kms~ 1 (/* = 0.09) and 
a C 2 = 579 km s -1 (/* = 0.02). In the M*-er C 2 plane 
(bottom panel of Fig. [T|) these models (represented by 
the left-most symbols in the diagram) lie well above 
the local observed correlation ((2]). We have little in- 
form ation on the veloci t y dis p ersion of high-z ETGs 
Csee iCenarro fc Truiillol 120091: iCappellari et all 120091 : 



Ivan Dokkum. Kriek. fc Franxl 120091 . hereafter ~vDKF09) : 
for 1255-0, an ETG at z = 2.186 with stellar mass 
M* ~ 3 x 1O U M (after conversion to Sa lpeter IMF) 
and h alf-light radius R c — 0.78 ± 0.17 kpc (Kric k~e~t al.l 
2009), vDKF09 tentatively measured projected velocity 



dispersion 510^g5 5 kms 1 , consistent within the uncer- 
tainties with cr e 2 of our seed galaxy models (see Fig. [TJ 
bottom panel). Given the uncertainties on <7 e 2, also the 
structural parameter c* is poorly constrained for high-z 
ETGs. Our seed-galaxy models — which were constructed 
to lie exactly on the lensing Mass Plane (NTB09) — have 
c* = 0.11 for /» = 0.09 and c* = 0.13 for /, = 0.02 
(left-most symbols in the top panel of Fig. [T]), consistent, 
for M r ~ 3 x 10 n M o , with the correlation © found 
for SLACS galaxies, and also with the highly uncertain 
value of c* estimated for the high-z ETG 1255-0. 

In Fig. [1] the points with log M*/Mq > 11.6 show the 
location of the merger remnants of all steps of the hi- 
erarchies in the M*-c*, M*-R e and M*-<j C 2 planes. For 
the adopted values of the length and mass scales, the 
timescale to complete one step of the merging hierar- 
chies is in the range 0.5 — 1.7 Gyr, so all hierarchies can 
be comfortably completed in the ~ 10 Gyr elapsing from 
z ~ 2 to the present. The location of the merger rem- 
nants with respect to the local scaling laws is quite de- 
pendent on the details of their merging history, but it 
is apparent that all the merger remnants lie significantly 
below the — R c relation and above the — er C 2 re- 
lation. On average, minor mergers increase more the 
galaxy size than major mergers, con s istent with the re- 
sults of lNaab. Johansson. &; Ostrikerl (|2009f ). but we also 
find that different minor-merging events produce rem- 
nants with significantly different R c and cr 2 for similar 
M*, even when the initial conditions differ just by the 
relative phase-space distribution of the centres of mass 
of the satellites and the central galaxy (see the right-most 
pairs of stars in Fig. [T]). 

The behavior of our minor-merger hierarchies in the 
M*-cr e 2 plane deserves some discussion. When a galaxy 
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grows by accreting satellites on almost parabolic orbits 
(as is the case for our minor-merger steps; see NTB09), 
the virial velocity dispersion of the remnant is lower 
tha n that of the progenitor if the mass loss is negligi- 
ble (iBezanson et al.1 120091 : iNaab. Johansson, fc Ostrikerl 
2009). The fact that in some of our minor-merger 
steps the remnant has slightly higher projected (and 
virial) velocity dispersion than the progenitor must be 
ascribed to the significant loss of luminous (see Fig. [T]) 
and dark (see NTB09) matter, which is known to 
have the effect of increasin g the velocity dispersion (see 
iNipoti. Londrillo. fc CiotHl2003l ). 

In general the merger remnants lie close to the FP- 
like relation M*-c* , consistent with previous studies con- 
sidering the FP (iGonzalez-Garcia fc van Albadal 120031: 
Nipoti. Londrillo. fc Ciottil 120031 : iBovlan-Kolchin et al I 
20061 ) or the lensing Mass Plane (NTB09). An excep- 
tion is the minor-merging hierarchy with /» = 0.02 seed 
galaxy (crosses in the diagrams), in which c* increases 
significantly with stellar mass, reflecting the relatively 
low M» and high cr e2 due to substantial mass loss. 

Altogether, considering both major and minor merger 
hierarchies, our results confirm that dry mergers bring 
compact ETGs closer to the observed correlations. How- 
ever, quantitatively, the process is not efficient enough 
to explain the observed size evolution. We note that the 
final merger remnants (the end-products of the last steps 
of the hierarchies) have exceptionally high stellar masses 
(M* >10 12 M Q ), so one cannot invoke more growth via 
dry mergers to obtain systems with surface mass den- 
sity as low as observed in local galaxies. Furthermore, 
velocity dispersion is approximately conserved along the 
merger hierarchies, and objects with <r C 2 ~ 500 kms -1 
are not observed in the local universe. 

4. CONSTRAINTS ON THE STELLAR-MASS ASSEMBLY 
HISTORY OF PRESENT-DAY EARLY-TYPE GALAXIES 

Taken at face value, the stellar masses and sizes of the 
high-z ETGs mean that these systems cannot evolve into 
present-day ETGs via dry mergers. However, possible 
observational biases in the high-z measur ements (for in- 
stance , that the sizes are underestimated; iHopkins et al.1 
2009bj) could make the discrepancy between high and 
low redshift systems smaller. Thus, if one only con- 
siders average properties, dry mergers may perhaps be 
still considered a viable evolutionary mechanism. How- 
ever, recovering the average size of local ETGs is not 
sufficient: additional constraints come from the small in- 
trinsic scatter of observed stellar-mass scaling relations, 
which constrains the dry-merger mass assembly history 
of present-day ETGs, independently of the actual com- 
pactness of their high-z counterparts. A first attempt 
of constraining the merging history of ETGs using the 
sc atter in the stellar -mass-size relation was performed 
by iShen et alj (|2003l ). who — on the basis of an analytic 
parametrization of the size growth due to dry mergers — 
found that a certain degree of fine tuning in the orbital 
parameters of the galaxy encounters is required in or- 
der to reproduce the local -M*-i? e relation of their SDSS 
sample of ETGs. We take here a step further, by taking 
advantage of N-body merger simulations and up-to-date 
stellar-mass scaling laws, including information on stellar 
velocity dispersion. 

We quantify the growth of galaxies by dry merging by 



approximating the dependence of the structural parame- 
ters on stellar mass with power-law functions (see NTB09 

for the lensing-mass analogs): R c oc M* R , cr o2 oc M* 

and c* oc M* c . To each step of a hierarchy we associate 
a value of Ar by calculating (3r = Alogi? G /AlogM», 
where A indicates variation between the remnant and 
the progenitor. The resulting 30 values of Ar are dis- 
tributed with average (/3b) — 1-09 and standard devia- 
tion 5(3r = 0.29. Doing the same for (3 a and (3 C , we get 
(Ar) = 0.07, 5(3„ = 0.11 and (A) = 0.22, 5(3 C = 0.13. 
Assuming that the galaxy merging history is character- 
ized by slopes (3 normally distributed with these averages 
and standard deviations, we can use these quantities to 
describe how dry mergers move galaxies in the planes 
M*-c*, M»-i? c and M^-a e 2- the relatively large values of 
5[3/(3 we find mean that dry mergers introduce substan- 
tial scatter in these planes. 

For the dry-merging scenario to be viable it must not 
introduce in the scaling relations more scatter than ob- 
served. We can use this condition to constrain the prop- 
erties of the candidate progenitors of local ETGs, as- 
suming as working hypothesis that they formed by dry 
merging. Let us consider, for example, local ETGs with 
given stellar mass in the range 10 11 — 10 12 M Q probed 
by the SLACS sample: in Fig.[2]the distributions (within 
l-cr) of c», i? c and <7 e 2 for these galaxies are represented 
by the vertical bars, while the loci in the M*-c„, M*-R e 
and M*-o~ e 2 planes of their allowed dry-merging progen- 
itors are shown by the shaded areas. Galaxies lying out 
of these areas, growing by major or minor dry mering, 
would populate regions of the parameter space inconsis- 
tent with the local observed distributions at the reference 
mass M*. It is apparent that the strongest constraint 
comes from the M»-i? plane (central panel), where the 
allowed region for progenitors spans ~ 0.26 dex in M*. 
This means that a growth in stellar mass by dry merg- 
ing by more than a factor of ~ 1.8 (and in size by more 
than a factor of ~ 1.9) is excluded on the basis of our 
set of simulations. In other words, present-day massive 
(M» ~ 10 11 - 10 12 M Q ) ETGs cannot have assembled, on 
average, more than ~ 45% of their stellar mass through 
dry merging. Quantitatively, the constraints on the mass 
assembly history depend on the specific choice of the or- 
bital parameters and mass ratios of the galaxies: in our 
set of merging hierarchy simulations there are more ma- 
jor than minor mergers, and we explore a wide range of 
orbital parameters. One might get different numbers by 
considering more minor-merging simulations or exclud- 
ing parts of the orbital parameter space. For instance, 
if we consider only the 8 minor merger steps, we get 
(As) = 1-30 and 8(3r = 0.40, which gives an upper limit 
of ~ 35% to the fraction of stellar mass assembled via 
dry mergers by local ETGs. In any case, it is clear that 
a remarkable degree of fine tuning is required to repro- 
duce the tightness of the local scaling relations with dry 
mergers. 

5. CONCLUSIONS 

Dry merging makes galaxies less compact. For this 
reason it has been proposed to be the key to reconcile 
the compactness of high-z ETGs with the lower stellar- 
mass surface-density of local ETGs. 

Our quantitative analysis based on N-body simulations 
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shows that it is very hard to reproduce with dry mergers 
the growth in size of a factor of ~ 5 between z > 2 and 
today reported in the literature for M* ~ 10 11 M Q ETGs, 
because the average size growth produced by dry merg- 
ing is not sufficient. The tightness of the local scaling 
laws of ETGs provides an additional and stringent con- 
straint for the dry-merging scenario. As we have shown, 
generic dry mergers tend to increase the scatter rather 
than decrease it because they cannot preserve simulta- 
neously the correlations between stellar mass, velocity 
dispersion and effective radius. On the basis of our set 
of major and minor merging simulations, we conclude 
that — independent of the actual compactness of higher- 
z ETGs — typical present-day massive ETGs did not as- 
semble more than ~ 45% of their stellar mass and grew 
more than a factor ~ 1.9 in size via dry merging. A 
similar upper limit is found also for the total (dark plus 
luminous) mass assembly on the basis of the lensing-mass 
scaling laws (NTB09). 

Is this then the end of dry mergers as a major evolu- 
tionary mechanism? The only way out could be provided 
by a high degree of fine tuning of the mix of progenitors 
and orbits. Although not a natural scenario in our opin- 
ion, this fine tuning could perhaps reproduce such tight 



scaling relations and thus save the day. A systematic 
exploration of cosmologically motivated merging hierar- 
chies is necessary to find out whether such an extreme 
fine tuning actually occurs. 

The question of what drives the apparent size evolution 
of ETGs remains open. A possibility is that a combina- 
tion of different astrophysic al processes and obser vational 
biases can do the job (e.g. iHopkins et aDl2009d) . but it 
is not excluded that the answer is to be found explor- 
ing more exotic scena rios, such as those assuming non- 
stand ard dark matter (|Lee 2008; Ferr er. Nipoti. fc Ettoril 
2009). 
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